STRONG SOLUTIONS OF STOCHASTIC EQUATIONS WITH 
RANK-BASED COEFFICIENTS 



TOMOYUKI ICHIBA, IOANNIS KARATZAS AND MYKHAYLO SHKOLNIKOV 

Abstract. We study finite and countably infinite systems of stochastic differential 
equations, in which the drift and diffusion coefficients of each component (particle) 
are determined by its rank in the vector of all components of the solution. We 
show that strong existence and uniqueness hold until the first time three particles 
collide. Motivated by this result, we improve significantly the existing conditions 
for the absence of such triple collisions in the case of finite-dimensional systems, 
and provide the first condition of this type for systems with a countable infinity of 
particles. 



1. Introduction 
We study the following system of stochastic differential equations: 

(1.1) dXi(t) = l{Xi®=x U) (t)} 8j dt + Yl 1 {X(t)=x U) (t)} <Tj dWi{t) 

for i G I . Here 7 = {l,...,n} for some n G N , or 7 = N ; 5j , j G 7 are real 
constants; <7j, j G 7 are strictly positive real constants; (Wj : i G 7) is a system of 
independent standard Brownian motions; and 

(1.2) X {1) (t) < X m (t) < X (3) (t) < . . . 

is the ordered particle configuration at time t . In addition, we let the initial config- 
uration be deterministic and satisfy 

(1.3) Xi(0) < X 2 (0) < X 3 (0) < ... . 

Ties in the ordered particle configuration are resolved in accordance with the initial 
ranking of particles; for instance, we set X^(t) = X{(t) , i = l,...,n whenever 
X 1 (t) = ... = X n (t) . We shall write X for (X { : i G 7) and W for (W t : i G 7). 

In the case 7 = N we work under the following assumption. 

Assumption 1. 7/7 = N, we assume that there is an integer MgN such that we 
have 

(1-4) 5 M = S M +i = ■ ■ ■ , 

(1.5) a M = a M +i = ■■■ ■ 



Date: September 20, 2011. 

2000 Mathematics Subject Classification. 60H10, 60J60, 60J65. 

Key words and phrases. Stochastic differential equations, strong existence, strong uniqueness, 
equations with rank-based coefficients, Brownian particles, triple collisions. 

This research was partially supported by NSF grants DMS-08-06211 and DMS-09-05754. 



2 



TOMOYUKI ICHIBA, IOANNIS KARATZAS AND MYKHAYLO SHKOLNIKOV 



Moreover, we assume that there exist constants 71 > 0, 72 G M such that 

(1.6) Xi{0) > 71^ + 72, i = 1,2,... . 

For the case I — {1, ... ,n}, the main result of [2] implies that the system (jl.ip has 
a weak solution ("weak existence"), which is unique in the sense of the probability 
distribution ("weak uniqueness"); the strict positivity of the diffusion coefficients is 
crucial here. In the case / = N , a slight modification of the proof of Proposition 
3.1 in [TB] shows that weak existence and weak uniqueness hold for the system (II. ip 
under the Assumption [TJ see Proposition H] below. 

In this paper, we investigate the questions of existence of a strong solution ( "strong 
existence" ) and of pathwise uniqueness ( "strong uniqueness" ) in both cases. Due to 
the discontinuity of the diffusion functions 

<Ti(x) = ^ 1 {xr=x U) } o-j, x E M n , % = 1, . . . , n 

in (11. ip . general results on strong existence and strong uniqueness, which rely on the 
regularity of the diffusion coefficients, do not apply even when / is finite. 

In order to construct a strong solution to the system (11.11) in the case / = 
{1, . . . ,n} , we rely heavily on the results of the recent article [8]; this paper deals 
with the case n = 2 and establishes strength and pathwise uniqueness for the solu- 
tion of the resulting system (11.11) (actually, even when one of the diffusion coefficients 
vanishes, but not both). The idea, then, is to put together paths of the strong solu- 
tions found in [5] for two particles, to obtain the strong solution in the case n > 2 of 
several particles. 

This is possible and the approach is viable, however, only when the particle system 
in (11.11) does not exhibit triple collisions, that is, when the event 

(1.7) {3 t > : Xi(t) = Xj{t) = X k {t) for some i < j < k } 

has zero probability for the state process X in the weak solution of the system (11.11) . 
We provide new, necessary and sufficient conditions for the absence of triple collisions 
in the case / = {1, . . . , n} ; and develop the first such conditions in the case I = N . 

To formulate our main results we shall need the following Definition [Tj as well as 
Conditions [T] and [2] below. 

Definition 1. A finite or infinite sequence (a\, 0,2, ■ ■ ■) is called concave, if for every 
three consecutive elements aj, aj + i, a i+2 we have 

1 / \ 
a i+ i > — (ai + a i+2 ) . 

Condition 1. The sequence (o~f, erf, • • •) is concave. 

Condition 2. Either I = {1, 2, . . . , n} and the sequence (0, af, erf , . . . , cr^,0) is 
concave; or I = N and the sequence (0, o~f, erf, . . .) is concave. 

Our main results read as follows. 
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Theorem 1. Consider the particle system in ( II. ip and, if I = N, let Assumption^] 
be satisfied. 

If the diffusion coefficients satisfy Condition HJ then the unique weak solution of 
( II. ip /ias no tnp/e collisions; that is, the event ( I1.7P /ias zero probability. On the 
other hand, if Condition^ fails, then the event (ll.7p /ias positive probability. 

Theorem 2. Consider the particle system in ( II. ip anrf, if I = N, let Assumption^] 
be satisfied. Introduce the first time of a triple collision, namely 

(1.8) r := inf{ * > | 3 i < j < k : X^t) = X^t) = X k (t)} . 

Then the system ( II. ip has a unique strong solution up to time r . 

In particular, if Condition^ is also satisfied, then there is a unique strong solution 
of the system (II. ip defined for all t > 0. 

Remark 1. Theorem [T] can be recast as saying that Condition [T] is necessary, and 
Condition [2] sufficient, for the absence of triple collisions. A condition that is both 
necessary and sufficient for the absence of triple collisions, has yet to be determined. 
So far, this question is completely resolved only in the case n = 3, in which the 
results of Varadhan & Williams [IT] imply that Condition [1] is both necessary 
and sufficient for the absence of triple collisions; see the proofs of Lemma [6] and 
Theorem [1] below for more details. 

The remaining gap between Condition [T] and Condition [2J is due to the following 
reason. By an inductive argument, we reduce the statement of Theorem [1] to the 
problem studied in De Blassie [7J. However, the (sharp) criterion given there 
involves the invariant distribution of the projection of a certain diffusion process in 
a Euclidean space on the unit sphere. Due to the lack of rotational symmetries in 
our situation, it is however not clear how to analyze this invariant distribution. For 
this reason, we simplify the condition in [7J to a checkable sufficient condition in 
Proposition \5\ below, sacrificing its sharpness at this point. 

Theorem [2] leaves open the questions of whether a strong solution continues to 
exist, and of whether pathwise uniqueness continues to hold, after a triple collision 
(we know from the work of Bass & Pardoux [2] that a weak solution exists after 
such triple collisions, and is unique in distribution). At the moment, we conjecture 
that strong solutions fail to exist beyond the time of the first triple collision, but this 
problem remains open and will have to be settled in future work. □ 

Questions regarding the presence or absence of triple or higher-order collisions in 
multidimensional diffusions have been addressed by several authors; in addition to 
the example of section 3 in Bass & Pardoux [2], one should consult the works 
by Friedman [9] (see also [10], chapter 11), Ramasubramanian [13], [H], De 
Blassie [7J, [6] and Ichiba & Karatzas [11] • In the papers [3] and [I], Cepa & 
Lepingle consider systems of Brownian particles with repulsive forces of electrostatic 
type, and show absence of triple collisions under conditions of sufficiently strong 
repulsion. 

Preview: The present paper is organized as follows. Section [2] collects a number of 
preliminaries, most notably results from [1], [12], [7J and [16] that are crucial in our 
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context. Section [3] deals with the absence of triple collisions under Conditional and 
with the proof of Theorem [TJ which represents a significant improvement over our 
earlier results in [TT]. In particular, we provide here new necessary conditions and 
new sufficient conditions for the absence of triple collisions in the case of a finite 
number of paticles, and develop the first such conditions for a countable infinity of 
particles. 

Section H] is devoted to the proof of Theorem [2J We start by setting up an inductive 
procedure, which "bootstraps" the strength of the solution to the system of equations 
01. ip that was established recently by Fernholz, Ichiba, Karatzas & Prokaj 
[8] for the case n = 2 of two particles - first to the case n = 3 of three particles; 
then to the case of an arbitrary, finite number n of particles; and finally, building 
on results of Shkolnikov [16], to the case of a countable infinity of particles. 



2. Preliminaries 



We start with some preliminaries on the weak solution of the system ( II. ip . when 
/ = {1, . . . , n}. First, we recall the dynamics of the ordered particles X^, . . . , X( n ) 
in the system (11. ip from section 3 in pQ and section 4 of [12]; once again, the strict 
positivity of the diffusion coefficients is crucial for these results. 

As in those papers, we shall denote by A J ' -1 ' J '(i), j = 2, . . . , n the local times (nor- 
malized according to Tanaka's formula) accumulated at the origin by the nonnegative 
semimar t ingales 



(2.1) Yj-iK-) '■= 

over the inteval [0,t] , and set A ' 1 ) 



X, 



X, 



A 



n,n+l 



(•) 



J ~- 

0. 



n 



Proposition 3. Set I = {1, . . . , n} and let (X, W) be a weak solution of the system 
(II. ip . Then there exist independent standard Brownian motions j3\, . . . , j3 n such that 

(2.2) dX U) (t) = 8 j dt + Gj d^(t) + - (dA J ' _lj '(f) - dA ij+1 (t)) , t > 
holds for all j € / . 

It was observed in [1] that Proposition [3] permits the identification of the process 

. . , Xua) as a Reflected Brownian Motion (RBM for short) 

{(xi,. . . ,x n ) E W 1 : x 1 < ...<x n } 



of ordered particles (Xm , 
in the wedge 

(2.3) QU:= 
with reflection matrix 



( -\ o 



\ 



i 

2 











1 

2 





That is, the process (^(i), • • • , X^) behaves as an n-dimensional standard Brownian 
motion in the interior of the wedge 2U , and is obliquely reflected on the faces {xi = 
Xi + i}, i = 1,. .. ,ra — l of QU. The directions of reflection are specified by the columns 
ti, i — 1, . . . , n — 1 of the matrix . 
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Occasionally, it will be more convenient to consider instead of the process of the 
ordered particles (X^, • • • , -X"(n)) the process of spacings (gaps) 

(2.4) Y := (X( 2 ) - . . . , X( n) - X( n -i)) 
as in ( 12. ip . From Proposition [31 we have the dynamics 

(2.5) dl^-i(t) = - df + o 3 dft(f) - a,-! dft_i(f) 



j (dA j ' j+1 (t) + dA^' -2 ^- 1 ^)) + dA J '- lj (t) , t > 



for the spacings of (12. ip with j = 2, . . . , n . Thus, the process Y is an RBM in the 
(n — l)-dimensional orthant (M+) Tl_1 with reflection matrix 



K :-- 



( 1 -\ \ 
-| 1 -I 

o -\ i ••. 



\ o o ■-. •-. J 

For a detailed summary of many results on Brownian motions with oblique reflection 
in the orthant, we refer to the excellent survey article [T8] . 

For further reference we make the following simple observation. The event in (I1.7P 
can be reformulated as 

(2.6) { 3 t > : Yi(t) = Y i+1 (t) = for some 1 < i < n - 1} , 

so that the presence or absence of triple collisions is an intrinsic property of the 
spacings process Y. 

Next, we let I = N and construct the weak solution to (11. ip along the lines of the 
proof of Proposition 3.1 in [IB] , under the Assumption [TJ 

Proposition 4. Let I = N and /et ^sstimpizonll] 6e satisfied. 

There exists then a weak solution (X,W) of the system (II. ip . which is unique in 
distribution. 

Moreover, after enlarging the probability space if necessary, we can find stopping 
times = kq < K\ < . . ., integers M = n(0) < n(l) < . . . and weaA; solutions X^ k \ 
k > of the system (11.11) with I = {1, . . . , n(fc)} swc/i i/iat 

(2.7) = X^\t), t e [« fc , K fc+1 ], 1 < z < n(A;), 

(2.8) X,(t) = Xi(0) + Sit + <TiWi(t), t e [K fc , K fc+ i], z > n(A;) + 1 

and, for each k > , the processes X^ and (W n ^)+i, W n ^)+2-, • • • ) ar & independent. 

Proof. The proof of Proposition 3.1 in [16] carries over mutatis mutandis to the 
situation here. We only need to replace the a priori estimate on the expected number 
of particles in an interval of the form (— oo, x] at a time t > by 



sup P ( Aj(0) - max |^| • t - sup / q(u) dW, 

~^ ?(•) V i 0<s<tj 



i(u) < X I < oo 



where the supremum is taken over all progressively measurable processes c(-) adapted 
to the filtration on the underlying probability space, which take values in the interval 
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[min, Gi , max, Oi ] . The finiteness of the series follows from the L 2 - version of Doob's 
maximal inequality for continuous martingales and the condition (jl.6p . □ 

Finally, as a preparation for the proof of Theorem [H we state a special case of the 
main result of De Blassie in [7J. 

Proposition 5. Let EI fee a finite- dimensional Euclidean space and a : EI — > EI 2 be 

a mapping satisfying 

a{x) = a f° r al1 {°} 

where || ■ || 2 denotes the Euclidean norm on EI, and suppose that the set of disconti- 
nuity points of er(-) on the unit sphere {x G H : ||x||2 = 1} has surface measure zero. 
In addition, let W be an H— valued standard Brownian motion, and suppose that the 
martingale problem corresponding to the H— valued stochastic differential equation 

(2.9) dV(t) = a(V(t)) dW(t) , 1/(0) G H \ {0} 
is well-posed. If the condition 

(2.10) inf ( 7 4^U > 2 

is satisfied, then we have 

¥(V(t) ^0, V t G [0, oo) ) = 1 . 

Here a(-) = <r(-)'cr(-) is the diffusion matrix of V, tr denotes the trace operator, 
and (• , •) is the Euclidean scalar product on EI . 

Proof. It suffices to note that (12.101) implies the condition 

inf B(x) > 1 

xeH 

INIl2=l 

in the notation of equation (1.9) in [7]. Thus, the result is a special case of Theorem 
1.1 (i) in [7J. □ 

3. Triple collisions 

The two main steps in the proof of Theorem [1] are provided by the following two 
lemmas. 

Lemma 6. Let I = {1, . . . ,n} with an integer n > 3 , and suppose that Condition 
[2] holds. Then the first time of a triple collision r, defined in (jl.8p , must satisfy 

(3.1) r = V 
with probability one, where 

(3.2) v := m£{t > : X 1 (t) = X 2 (t) = ... = X n (t)} . 



STRONG SOLUTIONS 



7 



Lemma 7. Let I = {1, . . . ,n} with an n > 3 and suppose that Condition^ holds. 
Then the first time of a triple collision r, defined in (jl.8p . must satisfy r = oo with 
probability one. 

We shall prove both lemmas simultaneously, by induction over n . 

Proof of Lemmas [6] and [TJ Step A. First, we note that for every T £ (0, oo) 
we can make a Girsanov change of measure such that the processes (Sit + o~i(3i(t), 
t £ [0,T]), i = 1, . . . ,n become independent standard Brownian motions under the 
new measure. Here, the Brownian motions . . . , /3 n are defined as in Proposition 

El 

Since almost-sure events retain this property under an absolutely continuous change 
of measure, and since T £ (0, oo) can be chosen arbitrarily, it suffices to prove the 
almost sure absence of triple collisions under the new measure. For notational sim- 
plicity, we prefer to assume S± — . . . — S n — from the start. Section 2.2 in [TT] can 
be consulted for a more detailed exposition of the same argument. 

Step B. We proceed with the inductive argument. For n = 3, we deduce from 
Proposition 3 of [TT] that 

(3.3) r = rj = oo 

holds with probability one, in the notation of (jl.8p and f)3.2p . In fact, this is a 
consequence of Theorem 2.2 in [17] for the reflected Brownian motion Y (see the 
proof of Proposition 3 in [TT] for more details). 

Step C. Now, fix an m > 4 and suppose that Lemmas [6] and [7J hold for all 3 < n < m. 
We will first show that Lemma |6] must hold for n = m as well. To this end, we define 
for each < s < 1 the stopping time 

(3.4) Ve := inf{t > : ||F(t)|| 2 < e or ||F(t)|| 2 > e" 1 }, 

where we have written || • || 2 for the usual Euclidean norm. 
We claim that, for all < e < 1 , the comparison 

t > rj e holds with probability one . 

If we can prove this claim, then we will be able to conclude that either r = rj = oo 
or t = 7] = lim £ ^o f] £ < oo must hold, and this will yield Lemma El 

To prove the claim, we deploy an argument similar to the one on pages 471-475 in 
[j~9] . As there, we consider the local behavior of the RBM of interest (in our case Y) 
on the compact sets 

(3.5) K £ := {ye(R + ) n - 1 : e < \\y\\ 2 < e' 1 } , < e < 1 . 

For each y £ (M + )™~ 1 \ {0} , we define an open set U(y) of (M + ) n ~ 1 such that 

(3.6) 3 1 < j < n - 1, 5 > : z 5 > S for all z £ U(y) , 
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and let j(y) be a number in {1, . . . , n — 1} as in (13 .6p . From the semimartingale 
decomposition of Y in (I2.5p . we see that if we start Y in the set U(y) for some 
y G (M + ) n_1 \ {0} , then we can write 

(Yj{t A Cum) ■ 1 < j < <y)) = (Y!(t A Cu {y) ) : 1 < J < a(j/)), 

(Y^t A Cu(y)) : b(y) < j < n - l) = A &,(„)) : 1 < j < c{y)) 

for all t > . Here we have set 

a(y) := j(y) - 1, b(y) := ){y) + 1, c(j/) := n - 1 - j (y) ; 

the process Y 1 is an RBM in the a,(y)— dimensional orthant; the process Y" is 
an RBM in the c(y)— dimensional orthant; and Cu(y) is the time that Y hits the 
boundary of U(y). In particular, the induction hypothesis implies 

( 3 - 7 ) r > Cu(y) , 

where y is such that Y(0) G U(y). 

Next, we fix an e G (0, 1) and cover the compact set K £ of (13. 5p by a finite number 
of open sets from the collection U(y) , ?/ G (M + ) n_1 \ {0} , say 

L 

(3.8) ^c[J U{y t ) . 

i=i 

Then we can find stopping times = Co < Ci < C2 < ■ • • of the form (u(y)j such that 
the path of the process Y(t A r) e ), t > can be decomposed into 

(3.9) Y(tA Ve ), Ck<t<Ck+u 



k > with the notation of (13.41) and with Y(t A r] E ) G U(yg) for some 1 < £ < L and 
all Cfc < t < Cfc+i- 

Using the strong Markov property of Y and the previous observation, one shows 
by induction over k that r > Ck A r/ e must hold with probability one, for all k > 0. 
By taking the limit fc — >• 00, we conclude that t > rj £ holds with probability one. 
Thus, outside of a set of probability zero, we must have 

(3.10) VeG (0,1) : r > rj e . 

This proves the claim and, thus, Lemma M for n = m. 

Step D. It remains to show that Lemma [7] holds for n = m. To this end, consider 
the centered process 

(n n \ 

X 1 (t) - n-^Xiit), ■ ■ ■ ,X n (t) - n-^Xiit) ) 
i=l i=l ' 

for < t < 00 . It is obvious that 

(3.12) 7] = m£{t > : V(t) = 0}. 

In addition, recalling that without loss of generality we have assumed 5% — ... — 
5 n = in Step A, we see from (II. ip that V is a diffusion process in the hyperplane 

(3.13) H = {xGM n :xi + ... + x n = 0} 
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with zero drift, and with diffusion function satisfying the conditions of Proposition 
Moreover, the martingale problem corresponding to the stochastic differential 
equation for the process V is well-posed, thanks to the main result of [2] and the 
non- degeneracy of the diffusion matrix of V on EI (which, in turn, follows from the 
strict positivity of the diffusion coefficients in ( II. ip ). Thus, we only need to check 
that condition (I2.10p is a consequence of Condition [21 since then Lemma [7] will follow 
from Proposition O 



To check (12.101) . we first compute the diagonal entries of the diffusion matrix a(-) 
of V in the coordinates of IR n to 

n 

(3.14) a*(l-2n- 1 )+n- 2 J2 a j > » = !,■■■,«, 

where the order depends on the ranking of the coordinates of V. Next, we note that 
the normal vector to EI is an eigenvector of «*(•) with eigenvalue . Thus, a(-) has 
an orthonormal eigenbasis over M n , which includes the unit normal vector to EI. It 
follows that the trace of the diffusion matrix cn(-) of F on EI coincides with the 
trace of the diffusion matrix of V in the coordinates of M. n , and is given by 

n n -, n 

Mo := £ (of (1 - an"') + n* £ oj) = 2zl . 

i=l j = l i=l 

Now, we estimate the denominator on the left-hand side of (I2.10p from above by 
the maximal eigenvalue of the diffusion matrix a(-) of V on EI. The latter is given 
by the maximal eigenvalue of the matrix 

(3.15) P'diag(a?,...,^)P, 

where P is the matrix representing the normal projection of vectors in M n onto EI , 
and diag(erf , . . . , a%) is the n x n diagonal matrix with diagonal entries af, . . . , <y\. 

Next, we observe that the spectral radius (and, thus, the maximal eigenvalue) of 
the matrix (I3.15P is given by 



(3.16) Mi := max 



a\x\ + ... + a 2 n x 2 n 

2 J- J_ t2 



zeK™\{0} V xi + . . . + X, 

xi+...+x n = 

This quantity can be further estimated from above as 

(3.17) M 2 := ma, ( ^ + ^ + 

x&."\{o} V xi + . . . + xi J 

x 1 +...+x n =0 v 1 n / 

where C is the maximal element of the set {erf, . . . , a 2 } and c is the second largest 
element of the same set. A careful optimization using Lagrange multipliers gives 

n ~ 1 ^ 1 
M 2 = C + — c , 

n n 

and shows that M 2 = Mi holds if and only if all elements in the set {erf, . . . , a%} are 
greater than or equal to c. 
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All in all, we conclude that the condition ( I2.10p . which proscribes triple collisions, 
amounts to 

(3.18) M Q >2M l 

and is satisfied, in particular, if the stronger inequality 



2 Mo 



1 n 11 

(3.19) M = — W > 2( — C+±c 

n ' V n n . 

i=i 

holds. We shall show (13.181) . or its stronger version (13.191) . by distinguishing the cases 
a\ C ^ o"^ and C = of (the case C = being completely analogous to the latter) . 

• In the first case (af ^ C ^ o^), the concavity of the sequence (0, of, . . . ,cr^,0) 
shows that its third-largest element is greater than or equal to \ C, whereas its 
fourth-largest element is greater than or equal to — C. Hence, the left-hand side of 
(I3.19P is greater than or equal to 

l 'c+c+±C+ic 



n V 2 3 

Plugging this into (I3.19P and simplifying, we see that it suffices to check (n — 3) c > 
-2^- C . The latter inequality, thus also ( 13 . 1 9[) and ( 12 . 1 0[) , holds for all n > 4 due to 
c > |C , a consequence of the concavity of the sequence (0, of, . . . , 0). 

• In the second case [C = af), we use the concavity of the sequence (0, af, . . . , <r^, 0) 
to deduce that its third-largest element is greater than or equal to -j C, whereas its 
fourth-largest element is greater than or equal to -r C . Hence, the left-hand side of 
( I3.19P is greater than or equal to 

[C + c+—C + — C) . 

n V 2 4 / 

Plugging this into ( 13 .191) . we conclude that it suffices to show 

Tl — 1 

(3.20) (n -3)c > ~^—C. 

Using c > 4- C (again, a consequence of the concavity of (0, af, . . . , cr^, 0) ), we 
observe: 

77 — 1 

(3.21) (n-3)c>^—C, Vii>4, 

with equality if and only if both c = -|- C and n = 4 hold. Moreover, the derivation 
of (I3.2ip shows that: either (13. 19[) holds (therefore, also (I3.18P ): or we have n = 
4, (of, of, erf, cr|) = (C, -| C, -j C, j C) and equality in ( 13 . 19|) . In this latter case, 
however, the inequality M 2 > M% in (13.171) is strict, so that (13. 18j) must hold. □ 

We can now combine our results, to prove Theorem [TJ 

Proof of Theorem [JJ Step 1. First, let I = {1, . . . ,n}. Then under Condition [2] 
there are no triple collisions, by virtue of Lemma [71 

Now, suppose that Condition [T] fails; that is, for some integer % — 2, . . . , n — 1 the 
comparison 

(3-22) a} - aU < ~ of 
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holds. Consider the weak solution X' = (X[, X 2 , X%) to the system (11. II) with / = 
{1, 2, 3} and the parameters 5i, Si + i, <Ti-\, <Ji, Oi+\ . Then, applying Theorem 2.2 
of [T7J as in the proof of Proposition 3 in [JTJ , we conclude that a triple collision of 
the particles X[, X 2 , Xg occurs with positive probability. It follows that there is a 
T G (0, oo) and a bounded, open subset U of the wedge {i 6 K 3 : X\ < %2 < 23} , 
such that the event 

{(x' {1) (t),x[ 2) (t),x[ 3) (t))eu, te[o,T}} 

n {3 i 6 [0,71: X[{t) = X' 2 {t) = X'S)} 

has positive probability for every initial condition in U. Along with the semimartin- 
gale decomposition of (12. 2p for the components of the process (X(i), . . . , X(„)) , this 
implies that the event 

(X _ 2) (t) ^ ^ X {i+2) (t), t G [0,2T] } 

n { (X (j „ 1) (t),X - ) (t),X (j+1) (t)) G U, t G [T,2T] } 

R { 3 t G [T, 2T] : *(,_!)(*) = X i3) (t) = X m) (t) } 

has positive probability. This completes the proof of Theorem [TJ for I = {1, . . . , n}. 

Step 2. Now, we turn to the case I — N and assume first that Condition [2] holds. 
We recall the notation in Proposition H] and observe that the event (11.71) is contained 
in the event 

(3.23) |J |J {X, n {t) = X i2 {t) = X i3 (t) for some te [« fc ,K fc+ i]}. 

Moreover, for every fixed %\ < i 2 < i% and k > , Proposition S] shows that there 
is a choice of n > M and a weak solution [x[ k \ • • • , Xn fc ^) of the system (11.11) with 
I — {1, ... ,n} and parameters 8x, . . . , 5 n , a±, . . . , a n , such that 

(3.24) Xi(t) = X^\t), t G [0, Kfc+i], 1 < i < n. 
Therefore, Step 1 of the present proof implies 

(3.25) P({ X h (t) = X i2 (t) = X h (t) for some t G [«jfc, }) = 

for any fixed i\ < i 2 < i 3 and k > 0. Thus, the event (11.71) has zero probability, as 
claimed. 

Next, suppose that Condition [TJ fails, and recall the definition of the constant M in 
Assumption [TJ In this case, there is an integer % — 2, . . . , M such that (I3.22j) holds. 
In particular, it follows from Step 1 of this proof that the weak solution X' of (II. ip 
with / = {1, 2, . . . , M + 1} and parameters Si, 5 2 , . . . , Sm+i, cr 2 , . . . , cr M+1 exhibits 
a triple collision with positive probability. In conjunction with the semimartingale 
decomposition of the process X', this shows that for every L > there is a T > 
such that the event 

{X[ M+1) (t) < X ( ' M+1) (0) — 27, T <t < 27} 
n {3 h < i 2 < i S) t G [T, 2T] : X^t) = X'^t) = X'Jt)} 
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has a positive probability. In addition, Proposition H] shows that, with probability 
one, the paths of the process are continuous and do not reach negative infinity 

in finite time. 

Putting these last two observations together, we see that there are exist real con- 
stants L > and T > , such that the event 

{ max(X {M+1) (t),X {M+1) (0) - L) < X (M+2) (t), t e [0,2T] } 

n { X {M+1) {t) < X (M+1) (0) - 2L, T < t < 2T } 

n { 3 i x < i 2 < i 3 < M + 1, t e [T,2T] : X h (t) = X i2 {t) = X i3 (t) } 

has a positive probability. In particular, the event (11.71) has a positive probability. □ 

4. Construction of strong solutions 

This section is devoted to the proof of Theorem [2J In the first subsection we 
explain our methodology in the special case / = {1,2,3}. The following subsection 
extends the construction of strong solutions to systems with any finite number of 
particles. Finally, in the last subsection we use the strong solutions in systems with 
finitely many particles to obtain the strong solution in the system with infinitely 
many particles. 

4.1. Systems with three particles. In this subsection we explain the construction 
of strong solutions in the case that there are only three particles. That is, we consider 
the system of SDEs 

3 3 

(4.1) dXi(t) = J2 l{x,(t)=x 0) (t)} Sj dt + J2 l{x l( i)=x (3)( i)} (Tj dWi(t) 

3=1 3=1 

with initial conditions satisfying ^i(O) < X 2 (0) < ^(O). As before, we define r to 
be the first time of a triple collision for a weak solution of the system (14.11) . To wit, 

(4.2) r = inf{t > : X 1 (t) = X 2 {t) = X 3 (t)}. 
In this setting our main result reads as follows. 
Proposition 8. Suppose 

2 2 \ 2 2 

O 2 -0- x > cr 3 - cr 2 . 

Then the system (14. ip admits a strong solution, which is pathwise unique. Moreover, 
if 

a 2 - a l < a 3 - 2 

holds, then a strong solution exists and is pathwise unique up to the triple collision 
time t. 

Proof. We will construct a strong solution to (14. ip by putting together paths of 
strong solutions to the system (11.11) with / = {1,2} and parameters Si, 5 2 , <Ti, a 2 and 
^2,S 3 ,a 2 ,a 3 , respectively. 

To this end, we introduce the following notation. For a closed time-interval [a, b] 
we denote by Z \- a < b ^> B > w \bi, b 2 ,ci,c 2 ) the M 2 -valued strong solution of the system (II. ip 
with drift parameters bi , b 2 and diffusion parameters ci , c 2 , driven by the indepen- 
dent standard Brownian motions B, W on the time-interval [a, b}. The latter strong 
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solution exists and is pathwise unique, on the strength of the results in section 5 of 
i- 

We can now define the stopping times = r < p < ti < P i < . . . and the desired 
strong solution on the intervals [r k , p k ], [pk,T~k+i], k>0 inductively by 

X M1) ([r k ,p k }) := {Z^ w ^' w ^(6 1 ,6 2 ,a 1 ,a 2 )) 1 , 

X*>Mn, Pk]) ■= {Z^ w ^ w ^(5 1: 5 2 ,a h a 2 ))^ , 
x n k (3)(t) := X^ (3) (r fc ) + 5 3 (t - r k ) + (T 3 W nk(3) (t) - (73^(3) (n), 

t e [r k ,p k ] 

p k := inf{t > r k : X nk(3) (t) = X nk(2) (t) or X nk(3) (t) = X nk(im } , 
Xe k (i)(t) := Xe k (i){p k ) + 5i(t - p k ) + <TiW 0k (i){t) - <TiW 0k (i)(pk) , 

t E [p k ,T k+1 ] 

Xe k{2) ([p k ,r k+1 }) := (Z^^ w ^ w ^(6 2 , 8,,^,^))^ 

X dk{3) ([ Pkl T k+1 }) := (Z^^W^W^^s.ffa.ffs))^ 

r k+1 := inf{t > p k : Xg k(2) (t) = X 6k{1) (t) or X 6k{3) (t) = Xg k[1) (t)}. 

For each k > , we have denoted here by ir k a permutation of the set {1,2,3} 
such that 

X nk(1) (T k ) < X^ k(2) (r k ) < X 7rfc(3) (r fc ) , 
and by 9 k a permutation of the set {1, 2, 3} such that 

Xe k {i)(Pk) < X ek{2) { Pk ) < X 9k ( 3) ( Pk ) . 

It is straightforward to check that the just constructed processes X\,X 2 ,X 3 are 
well-defined and form a strong solution to the system (14.11) up to the time 

(4.3) t := lim r k = lim p k . 

k— >oo k— >oo 

On the other hand, since the paths of the ranked processes X^, X^, X( 3 ) are 
continuous, we have 

X (1) (r ) = lim X (1) (r fc ) = lim X {2) {r k ) = X {2) (r). 

k — y oo k~ — yoo 

Moreover, an analogous computation yields X^(t) = X( 3 )(r ) . Thus, r > r. This 
proves the existence results of Proposition [HJ 

We now turn to the pathwise uniqueness of the solution. In the case a\ — o\ > 
a i ~ °\ i ^ follows from Theorem [T] that r = oo with probabilty one. In this case, 
pathwise uniqueness of the solution constructed above is a consequence of Theorem 
3.2 in [5]. The latter states that, for a system of stochastic differential equations with 
time-independent coefficients, strong existence in the presence of weak uniqueness 
implies pathwise uniqueness. 

In the case o\—a\ < a\—a\ , we recall from Theorem[T]that we have P(r < oo) > 0. 
In this case, we let (X(t),W(t)), < t < r and (X(t),W(t)), < t < ? be two 
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strong solutions of the equation (14.11) . where r is the first time of a triple collision in 
the particle system X = (Xi, X 2 , X3). By enlarging the probability space if necessary, 
we can extend X and X to weak solutions of the equation (14. ip . defined on the whole 
time interval [0, 00). Then, an application of Theorem 3.1 in [5] shows that the joint 
laws of the triples (r, X([0, r]), W([0, r])) and (f, X([0, ? ]), W([0, r ])) are the same. 
We can now proceed as in the proof of Theorem 3.2 in [S] to deduce r = r and 
X([0,t]) = X([0,t]) with probability one. □ 

4.2. Systems with finitely many particles. We now turn to the proof of Theorem 
[2] with I — {1, ... ,n}, where n > 3. Although the main idea behind the construction 
of the strong solution is the same as for n — 3, the proof is more involved here 
due to a more complicated pattern of collisions. For example, even in the absence 
of triple collisions, it is still possible to have collisions of the form Xi(t) = X 2 (t) , 
X 3 {t)=X 4 {t). 

Proof of Theorem [5] for I = {1, . . . , n}: Step 1. As in the proof of Proposition |H1 
we start by constructing a strong solution to the system (11.11) in an inductive manner. 
However, this time several layers of inductive constructions will be necessary. 

First, we recall the notation Z^ a > b ^> B,w (bi, b 2 , ci, c 2 ) for the strong solution of the 
system (II .ip with / = {1,2} and parameters 61,62,01,02, which is driven by the 
independent standard Brownian motions B, W on the time interval [a, 6]. The latter 
exists by the results of section 5 in [5]. Next, we define a sequence of stopping times 
= To < 7i < . . . , subsets A , A\, A 2 , ... of the set N with A = 0, and the desired 
strong solution of (II .ip on [r , t±], [r±, r 2 ], . . . inductively by 

j : {j-l,j}n4 = 0, te [r k ,r k+1 ], 

X nkU) [r kl r k+1 ] = ( Z ^M k u^ k u + v( 5j , 5 j+h * jt a j+1 )) v 
X Mj+1) [r k , r k+1 ] = (Z^^W^w^ y j} 5j+u aj} 

J e A k , 

r k+1 = inf{t > r k : X {j) (t) = X (j+1) (t) for some j <£ A k }, 

A k+ i = {1 < j < n - 1 : X {j) (r k+1 ) = (r k+1 ) } , 

where for each k > 0, ir k is a permutation of the set {1, . . . , n} for which 

x ir k (i){i~k) < ^7r fe (2)(u0 < • • • < X nk{n) (r k ). 

As in the case of n = 3 it is not hard to see that this defines a strong solution of the 
system (11. 1]) up to the time 

r' 1 ' : = lim r k . 

k— >oo 

We note at this point that, for n > 3, r' 1 ' = r is not necessarily true. 

Step 2. To proceed, for 1 < j < n — 1 we let (pj,k)k>i be a (possibly empty) 
subsequence of the sequence (r k ) k >i, which contains all the elements of the sequence 
(j k ) k >i for which X^{r k ) = X^ +1 )(r k ). Since for every k > there are at least two 



STRONG SOLUTIONS 



15 



sequences of the form (pj,k)k>i, which contain at least one element of {rk , Tk +i}, at 
least two of the sequences (pj,k)k>i are infinite. Thus, due to the continuity of the 
paths of the ordered particles, there exist 1 < ji < j'2 < n — 1 such that 

*Ch)(T M ) = % 2 )(r [1] ) = X (j2+1) (rW) 

holds. 

If J2 = ji + 1? we have a triple collision of X^, .X^'i+i) and X(j 1+2 ) at time rW. 
In this case, r < r^, and the existence results of the theorem readily follow. 

If h + 1 ) then we proceed with the construction of the processes X\, . . . , X n as 
in Step 1, but now starting at t = instead of t = 0. This gives us a new sequence 
of stopping times = Tq 1 ' < rf 1 < . . . defined analogously to the stopping times 
To < ti < . . . in Step 1. Next, we set = lim^oo and observe that we have 
constructed a strong solution to the system (II. ip up to the time r™. 

Again, either there is a triple collision at and the proof is complete, or we 
proceed with the construction to extend the solution up to a time r' 3 l Proceeding in 
the same manner and assuming that we do not encounter a triple collision, we end 
up with a strong solution up to time r'°°] := Imv^ooT^ . 

Now, assuming that a triple collision has not occured, we conclude that, for each 
fc£N, there exist 1 < ji(k) < ]2{k) < n — 1 such that ]2{k) > ji{k) + 1 and 

X (Mk))( T[k] ) = X (j2(k)){.T [k] ) = X {j2{k)+1) (T [k] ). 

Arguing as before, we conclude that there exist 1 < ji(oo) < ^2(00) < ^3(00) < 
.74(00) < n — 1 such that 

X(Moo))(r [oo] ) = % H+ i)(r H ), I = 1,2,3,4. 

Again, either there is a triple collision at time r^ 00 ', or we proceed with the con- 
struction of the strong solution until a time f at which 

X Cd f )= X Cn^ * = -2*, 

where K is such that 2 K > n/2 (unless a triple collision occurs before the time f). 
Now, the inequality 2 K > n/2 implies that at time f there must be a triple collision, 
which in turn implies r < f. Thus, we have constructed a strong solution up to the 
time r, as desired. 

Finally, pathwise uniqueness up to time r can be shown by the same arguments 
as in the case n = 3 (see the proof of Proposition [S] for the details). □ 

4.3. Systems with infinitely many particles. We can now use the results of 
the previous subsection to construct the unique strong solution for the system with 
infinitely many particles. 



Proof of Theorem [2] for I = N: Consider a probability space on which a system 
W = (Wi : % £ I) of independent standard Brownian motions is defined. Then the 
result of the previous subsection shows that for every n e N, there is a strong solution 
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to the system (II. ip with I = {1, . . . , n} and the parameters S\, . . . , 5 n , a±, . . . , a n 
which is defined up to the time 

(4.4) inf{t > : xf l \t) = xf\t) = X^\t) for some 1 < i < j < k < n}. 

In particular, it is defined for all t > if Condition [2] is satisfied. Now, we can 
construct a strong solution X of the system (II. ip up to time r by following the 
proof of Proposition 3.1 in [IB] , and making the following two modifications. We use 
the just described strong solutions X^ instead of the weak solutions used in [TB| 
and replace the a priori estimate for the expected number of particles on intervals 
of the form (— oo,x] at a time t > in the proof of Proposition 3.1 in [16] by the 
corresponding a priori estimate in the proof of our Proposition HI This proves strong 
existence. 

To prove strong uniqueness, let X' = (X- : i G /) be another strong solution 
of the system (11.11) . defined on the same probability space as X and adapted to 
the same Brownian filtration, and let r' be the corresponding first time of a triple 
collision. Moreover, we define the stopping times = k' < 4 < . . . and the sets 
{1, . . . , M} = I' C I[ C . . . inductively by 

4 +1 = (infO > 4 | 3 i G I f k ,j i I' k : X>(t) = X>(t)}) A r', 

I' k+1 = { % G / |3 < t < K ' k+1 : X[{t) = X'jit) for some j G l' k } 

and let = Kq < K\ < . . . and {1, . . . , M} — I C I\ C . . . be the corresponding 
quantities for the strong solution X. Then, by the strong uniqueness for the finite 
system, for each k G N the process 

A <),...,*[,,,(* A <)), t>0 

must be the unique strong solution of (11.11) with I = I' k and parameters Si, ... , S\p\, 
01, . . . , <7|j/| driven by W±, . . . , W\p\, stopped at n' k . Moreover, we have 

X[{t A 4) = X[{Q) + 5i(t A 4) + (TiWiit A 4) , t > 

for all i (fi I' k . The same arguments with X' replaced by X and induction over k give 

(4.5) 4 = K fe , I' k = I k , X(t A Kk ) = X'{t A 4) 

for all t > and k G N. Thus, by taking the limit k — > oo we obtain r = r' and 
I(tAr) =X'(t At'), t > 0. 

Note that the almost sure identities 

t = lim Kk and r' = lim 4 

fc— >oo fc— >oo 

follow from the fact that, with probability one, for every cc G K and t > there are 
finitely many particles on the interval (— oo,x] at time t in X and X' (due to the a 
priori estimate in the proof of Proposition H] and the Borel-Cantelli Lemma). □ 
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